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Abstract 

We  extend  the  linear  “stick-slip”  models  of  Doi-Edwards  and  Johnson-Stacer  to 
nonlinear  tube  reptation  models.  We  then  show  that  such  models,  when  combined  with 
probabilistic  formulations  allowing  distributions  of  relaxation  times,  provide  a  good 
description  of  dynamic  experiments  with  highly  filled  rubber  in  tensile  deformations. 

A  connection  to  other  applications  including  dielectric  polarization  and  reptation  in 
other  viscoelastic  materials  (e.g.,  living  tissue)  is  noted. 

Keywords:  viscoelastic,  hysteresis,  polymers,  molecular  models,  relaxation  times,  probability 
distributions,  uncertainty 


1  Introduction 

This  note  is  prompted  by  several  thrusts  in  our  research  efforts.  The  first  is  to  extend 
linear  reptation  models  for  polymeric  materials  to  models  incorporating  nonlinearities  and 
to  use  the  resulting  systems  to  explain  molecular  based  hysteresis  (e.g.,  via  internal  variable 
formulations).  A  second  direction  involves  exploration  of  multiscale  aspects  of  polymeric 
structural  modeling  with  uncertainty  at  the  molecular  (micro)  level.  We  do  this  in  the 
context  of  a  probabilistic  formulation  of  the  models  to  produce  a  suitable  overall  system 
(macro)  response  to  deformations.  The  ideas  are  illustrated  in  a  specific  application  to 
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Figure  1:  Experimental  stress-strain  curves  for  (1)  unfilled,  (2)  lightly  filled  and  (3)  highly 
filled  rubber  in  tensile  deformations 

highly  filled  rubber  exhibiting  significant  hysteresis  as  well  as  nonlinearity  in  tensile  and 
shear  deformations  as  depicted  in  Figure  1. 

Indeed,  the  ideas  reported  here  were  motivated  by  earlier  efforts  (summarized  in  [7] 
and  the  references  therein)  using  a  phenomenological  approach  to  filled  rubber  modeling 
where  we  encountered  significant  hysteresis.  The  phenomenological  approach  entailed  the 
use  of  Boltzmann  hysteresis  operators  with  nonlinearities  that  were  necessary  to  describe 
experimental  data.  Our  desire  to  understand  the  models  at  a  more  mechanistic  level  led  to 
the  linear  molecular  models  of  Doi  and  Edwards  [15]  and  Johnson  and  Stacer  [20].  Linear 
versions  of  those  formulations  did  not  lead  (see  [6])  to  the  types  of  nonlinear  hysteresis 
formulations  sought  after. 

The  focus  here  is  on  viscoelastic  polymers,  but  the  underlying  ideas  are  much  more  widely 
applicable  to  problems  in  biology  (living  tissue,  disease  pathogenesis  [5]),  dielectric  materials 
(polarization  effects),  industrial  fluids  (polymeric  melts),  and  ecological  migrations  (hidden 
or  internal  episodic  behaviors).  For  example,  treating  shear  waves  in  living  tissue  requires 
nonlinear  constitutive  laws  that  are  hysteretic  in  nature  ([3,  8,  17]).  Molecular  level  (internal 
strain)  formulations  also  utilize  multiple  relaxation  time  constants  precisely  such  as  those  in 
the  models  developed  below.  Multiple  relaxation  times  play  an  important  role  in  molecular 
based  (tube  reptation)  nonlinear  constitutive  models  for  the  flow  of  polyethylene  melts  ([11, 
12,  13,  18,  22]).  For  dielectric  materials,  it  has  been  known  for  some  time  ([14,  25,  26])  that 
one  needs  distributed  relaxation  times  in  polarization  models  for  heterogeneous  materials. 
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For  example,  it  is  desirable  to  represent  the  complex  dielectric  permittivity  ([14])  in  terms 
of  a  probability  density  (j)  over  relaxation  times  as  in 


£(UJ)  ~  £ca  +  (^s  —  £oo) 

This  can  be  written  in  a  more  familiar  form  as 


0(r)dr 
1  +  jur 


£(cj)  —  £oo  +  (^s  —  £oo) 


(j>(r)dT 

1  +  c o2t2 


j(£s  ^oo) 


a;r0(r)dr 
1  +  uj2t 2 


The  latter  formulations  are  analogous  to  proposed  forms  of  the  “elastic”  moduli  based  on 
experimental  observations  in  the  works  of  Ferry,  Andrews,  Ter  Haar  and  others  [2,  16,  23, 
24,  27].  We  describe  in  detail  the  incorporation  of  multiple  relaxation  times  in  our  models 
in  Section  4. 

In  Section  2  we  give  a  brief  review  of  the  linear  reptation  models  developed  by  Johnson 
and  Stacer  and  then  in  the  next  section  provide  the  details  on  our  extension  to  nonlinear 
versions  of  these  models.  In  Sections  4  and  5  we  explain  how  these  models  have  led  to  a 
molecular  based  fit  of  dynamic  models  to  experimental  data  for  highly  filled  rubber  in  tensile 
and  shear  deformations. 


2  Linear  reptation  models 


Tube  reptation  models  for  deformations  of  viscoelastic  polymers  were  introduced  by  Doi  and 
Edwards  [15],  and  were  further  developed  by  Johnson,  et.ah,  in  several  papers  [19,  20].  In 
this  section  we  give  a  brief  overview  of  their  models  and  assumptions  in  the  case  of  tensile 
deformations  since  they  serve  as  a  starting  point  for  the  derivation  of  our  nonlinear  reptation 
model  in  Section  3. 

The  Doi-Edwards  “stick-slip”  model  assumes  that  the  polymer  is  composed  of  chemically 
cross-linked  (CC)  tubes  that  contain  physically  constrained  (PC)  molecules.  At  the  time  an 
instantaneous  tensile  step-strain  is  applied,  the  entrapped  PC  molecules  stick  to  the  tube 
and  elongate  with  it,  but  then  they  contract  and  slip  back  close  to  their  original  length.  This 
provides  the  viscoelastic  character  of  the  material.  In  particular,  the  relaxation  is  modelled 
in  the  following  way.  Let  L(t)  denote  the  length  of  the  chemically  cross-linked  molecule, 
while  £{t)  stands  for  the  length  of  the  physically  constrained  one.  Assume  that  the  PC 
molecule  elongates  to  length  £*  due  to  the  applied  step-strain.  Then 

£(t)  =  £(0)  +  {£*  -  £(0))e~t/T  0  =  to<t<t1,  (2.1) 


where  r  is  the  relaxation  time  for  the  “slip”  motion,  and  t\  denotes  the  time  the  next  step- 
strain  is  applied.  It  is  assumed  in  [15]  and  in  the  continuum  realization  of  this  model  in  [20] 
that  the  “stick”  phase  of  the  motion  at  the  time  of  the  step-strain  deforms  the  PC  molecule 
proportionally  with  the  CC  tube  deformation,  i.e., 


A£i 

A  L, 


i  =  0, 1, ... . 


(2.2) 
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Here  Llt  and  £t  denote  the  length  of  the  CC  and  PC  molecules  at  time  U,  respectively,  while 
A L{  and  A C  stand  for  the  instantaneous  stretches,  i.e.,  L,  =  Lj_i  +  ALj_i,  £t  =  C-i  + AC-i- 
Using  this  assumption  we  can  write  (2.1)  as 


£(t)  =  £q  +  -^-A L0e  <//r,  0  —  t0  <  t  <  t\.  (2.3) 

Co 

Continuous  motion  is  approximated  by  a  series  of  step-strain  deformations  applied  in  At 
intervals,  where  At  — >  0.  This  leads  to  the  relaxation  equation 


or,  in  differential  form, 


Jo  L(s)  as 


di  _  4  _  (I  _  1  dL\ 
dt  t  \t  L  dt  J 


(2.4) 


When  this  model  is  assembled  into  a  3-D  continuum  model  of  a  solid,  one  defines  a  unit 
box  or  cell  at  each  material  point  oriented  by  the  principal  stretches  (Aic,  A2c,  A3C)  of  the  CC 
system  with  an  inside  box  with  parallel  sides  (Alp,  A2p,  A3p)  for  the  PC  system.  Stresses  are 
calculated  by  determining  how  the  strain  energy  function  W,  which  has  contribution  i'Ucc 
from  the  CC  box  and  Wpc  from  the  PC  system,  changes  with  respect  to  the  applied  stretches 
or  displacements  of  the  CC  system.  In  turn,  the  stretches  of  the  PC  box  are  treated  as 
internal  variables  depending  on  the  stretches  of  the  CC  system.  We  note  that  as  an  analogy 
to  (2.2)  we  have 


r)\  \ 

u/xjp  Jp  x 
dXjr  Xir 


(2.5) 


where  Sji  is  the  Kronecker  delta  function.  Thus  we  have  that  the  strain  energy  density 
function  is 


If  —  Wpc  +  IUCC  —  Wpc  (Aip,  A  2p,  A3p)  +  Wcc(Aic,  A2c,  A3c)  , 

and  the  internal  dynamics  (2.4)  (with  d0  —  1)  yields 

dXjp  _  1  _  (l _ ^dXic\ 

dt  t  \t  XiC  dt  J  ipi 


(2.6) 


(2.7) 


since  the  principal  directions  of  the  PC  and  CC  system  are  aligned.  The  Cauchy  stress  in 
the  principal  direction  ej  (where  e3  is  a  unit  vector  in  the  x3  direction)  is  given  by 


dW 


Tj  XjcdX 


-  P, 


JC 


(2.8) 
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where  P  is  the  hydrostatic  stress.  We  combine  (2.6)  and  (2.8)  with  the  assumption  that  A iP 
depend  on  A jC  as  internal  variables  and  obtain 


>  dWcc  *dWpcd  Xip 

Tj  jc  d\jc  +  irf-  d\ip  d\jc 

_  ^  dWcc  ^  dWpc  _  p 

nr.  r\  \  I  /v  in  /  \  x  JP 


JC  dXjc 


JP 


dX 


1 jp 


by  applying  (2.5).  If  we  choose  j  —  1  for  the  direction  of  loading,  we  have 


t2  =  t3  =  0  =  A 


d  Wn 


2  c 


dX 


+  A 


dW. , 


2c 


2  P 


pc 


dX 


P 


2  P 


(2.9) 


because  stresses  are  zero  on  the  sides  of  a  tensile  sample.  Since  P  can  be  determined  from 
(2.9),  we  obtain  that  the  tensile  Cauchy  stress  is 


T\  —  X 


'0  Wn 


1 C 


-A 


dWn 


\  dXic  dX2c 

while  the  engineering  stress  is  given  by 


+  A 


i  p 


'dwpc  dWp; 

—  a2 


dX 


i  p 


'2p 


dX 


2  P 


O'  1  = 


Tl 

Aic 


A 


2c 


dWa 

dX 


+ 


Aip  (dW, 


pc 


2c 


X 


1  c 


dX 


i  p 


x  owp; 


This  model  was  analyzed  in  [6]  for  incompressible  rubber  materials  undergoing  large  dynamic 
tensile  strains  with  a  particular  strain  energy  function  provided  in  [21]. 

Let  up,uc  denote  the  deformation  of  the  PC  and  CC  boxes,  respectively.  With  Aip  = 
1  +  dxup  and  Aic  =  1  +  dxuCl  (2.7)  implies  that 


d_ 

Ft 


(it  0xup ) 


i 

r 


(1  T  dxuc )  1 


d2uc 

dtdx 


(IT  dxup) . 


This  in  turn  yields  the  internal  dynamics 

1 

£\  H - E\  =  T  £i)/(1  T  s) 

T 

in  terms  of  the  infinitesimal  strains  ey  =  dxup  for  the  PC  system  and  £  =  dxuc  for  the 
CC  system.  It  is  shown  in  [6]  that  the  above  derivation  leads  to  an  approximation  of  the 
engineering  stress  in  the  form 

<7i(e,  ei)  w  (1248  -  1014^)^  +  1014$^  =  (1248  -  1014  ei)e  +  1014ei. 

ox  ox  ox 

This  can  be  combined  with  the  basic  model  for  the  longitudinal  vibrations  of  a  rubber  rod 
as  given,  for  example  in  [10],  by 

()~  /  j  cPii  c) 

]  =  F’  (210) 

where  F  is  the  applied  external  force,  Ac  is  the  cross-sectional  area  and  p  is  the  mass  density 
of  the  material. 
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3  Nonlinear  reptation  model 


In  this  section  we  present  a  nonlinear  extension  of  the  “stick-slip”  model  of  Doi  and  Edwards 
for  tensile  deformations.  A  crucial  assumption  throughout  the  derivation  in  the  previous 
section  was  that  the  elongation  of  the  PC  molecules  is  proportional  to  that  of  the  CC 
molecules  during  the  “stick”  phase  of  the  motion,  i.e.,  ^  or  5yi  as  in  (2.2) 

and  (2.5).  It  might  be  expected  that  for  certain  materials  and  large  strains  this  relationship 
is  not  strictly  linear,  but  rather,  it  is  described  by  a  general  nonlinear  function  /  that  may 
also  depend  on  time,  i.e., 


—  =  fit  -)  or  =  fit  —)6  - 

A  L  n,L)f  d\r 


(3.1) 


Thus  we  have  that  the  relaxation  equation  after  an  instantaneous  step-strain,  in  contrast 
with  (2.3),  is  given  by 

m  =  4  +  /(to,  ^)A L0e-(t-*o)/T,  0  =  to  <  t  <  h- 
Similarly,  for  t  e  (tm-i,tm),  m  >  1  we  have 

m— 1  n 


In  the  limit,  as  At  =  tr 


Thus 


m  =  t0  +  £  f(U,  f)AUe-^/\ 

i= 0 

tm _ i  tends  to  zero  we  obtain 


dt  =  £(t)  dL 

dt  ~n  ,L(t)’dt 


describes  the  continuous  motion  of  the  CC-PC  system  in  differential  form.  Now  we  take 
£o  =  l,  £  =  1  +  £i,  L  =  l  +  £  and  we  obtain  the  internal  dynamics 


£i  H — £i  —  £/(t,  y— — , 
r  1  +  e 


(3.2) 


As  in  the  previous  case  we  can  add  this  equation  to  the  general  model  for  the  longitudinal 
vibrations  of  a  rubber  rod  (2.10). 

Before  turning  to  the  engineering  stress  oq  and  describing  how  it  is  affected  by  the 
nonlinear  assumption  (3.1),  we  consider  approximations  of  (3.2).  Assuming  that  /  is  time 
invariant  and  expanding  it  in  a  Taylor  series,  we  obtain 


/( 


1  +  S\  , 

l  +  £  ' 


=  /(l) +  /'(!)( 

=  /(!)  +  /'(!) 


1  +  £\ 
l+£ 
£1  —  £ 
l  +  £ 


1  + 


/"(l) 


1  +  £\ 
1+  £ 


+  ... 


+ 


/"(!)  (£_ l  ~  g 
1  +  £ 


+  .  .  . 


=  70  +  fei  +  7i£  +  72e2  +  h.o.t.(ek,  k  >  3)  +  h.o.t.(£k,  k  >  2) 
+  h.o.t.(£k£{,  k,j  >  1). 


6 


If  ci  <<  e,  we  find  the  approximate  equation  for  (3.2)  given  by 

ii  +  -£1  =  e(7o  +  7ie  +  72£2)- 
T 


(3.3) 


We  note  that  this  approximation  corresponds  to  our  earlier  phenomenological  internal  vari¬ 
able  formulation  in  [10] 

•  ,1  d  /  \ 

£i  H — £i  —  xtIM6); 
t  at 

where  gv  is  a  cubic  polynomial. 

Now  we  turn  to  the  derivation  of  an  approximation  to  the  engineering  stress  ay  in  presence 
of  the  nonlinear  assumption  (3.1).  Let  us  make  similar  assumptions  to  those  that  we  made 
for  the  linear  model  in  Section  2.  Namely,  we  have  that  the  strain  energy  function  W  is 
given  as 


II  —  Wpc  +  Wcc  —  Wpc(Alp,  X 2P,  A3 ip)  +  Wcc( Aic,  A2c,  A 3c) ,  (3-4) 


so  that 


dW 


Tj  ~  XjcdX 


P, 


JC 


where  the  Aip  and  A,;c  are  the  principal  stretches  of  the  PC  and  CC  systems,  respectively. 
Additionally,  we  assume  that  /  in  (3.1)  is  independent  of  time,  i.e.,  fit, 6)  =  f{6).  Since  the 
A jP  depend  on  the  Aic  as  internal  variables,  the  Cauchy  stress  can  again  be  given  in  the  form 


T> 


=  A 


d  Wcc 
'JC  d\jc 


*dWpcd  xip 
jchdKdxjc  • 


(3.5) 


At  this  point  our  derivation  differs  from  the  one  in  the  previous  section.  By  (3.1) 


*9  A  ip  ..  /  Xjp .  ~ 

ox 7c  =  nxiJ  ij ’ 


so  (3.5)  yields 


-  >  SWcc  ,  ,  ±dWVc,,Kr^  „ 

Tj  —  ^jc  +-^ic2^  -  /(  ,  )$ij  P- 

OAjc  ^=1  (JAip  Ajc 


(3.6) 


If  j  =  1  denotes  the  direction  of  loading,  then  we  have 


n  x  dWcc  ,  x  dWpCif a2pn  „ 

r2  =  r3  =  0  =  A2c— —  +  A2c— —  /(— )  -  P 
oA2c  a2c 


(3.7) 


We  can  express  P  from  (3.7)  and  substitute  into  (3.6)  with  j  —  1  to  find 


,  d Wcc  ,  ,  dWpCftXlp, 

ti  =  Alc— - h  Aic- 


<9A 


lc 


<9  Aip  Aicy 


,  <9 Wcc  ,  .  «9WpC/.A2p 

A2c  VTT - H  A2cvp; — J  W — 

<9A2r  <9A2p  A 


2c 
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(3.8) 


Now  the  engineering  stress  a i  is  given  by 

=  dWcc  _  X tcdWcc  dWpcfA ip,  _  X2cdWpc  X2p 

dXi c  Aic  dX2c  dX\ p  Aic  Aic  dX2p  X2c 

The  first  two  terms  on  the  right  side  above  constitute  the  contribution  of  Wcc  to  the  engi¬ 
neering  stress  <7i,  while  the  last  two  terms  provide  the  contribution  of  Wpc.  We  assume  that 
the  energy  density  function  W  is  given  as  in  (3.4)  with 

Wcc  =  C\(I\  —  3)  +  C2(I2  —  3), 
wpc  =  c3(h  -  3)  +  C4(Il  -  3)2  +  C5(1 2  -  3)3, 

where  specific  values  of  the  coefficients  C\, . . . ,  C$  can  be  chosen  as  in  [6]  given  the  strain 
energy  function  suggested  by  Johnson  and  Stacer  based  on  experimental  data  [21],  Here  the 
strain  invariants  are 


h 


X2 

me 


+  A2C  +  A 


2 

3c’ 


X2  X2 

/'le/'2c 


+ 


+  ^2c^ 


2 

3c’ 


while 


h  —  X\p  +  A2„  +  A2 


'2  p 


3  p’ 


^lp^2p  +  ^lp-^3p 


+  ^2 


2 

3  P' 


We  also  impose  the  incompressibility  condition  in  the  principle  stretches  for  the  PC  and  CC 
systems,  respectively,  i.e.,  AipA2PA3P  =  AicA2CA3C  =  1.  We  can  compute  the  contribution  of 
Wcc  to  the  engineering  stress  exactly  the  same  way  as  in  [6] 


d  Wcc  X2cd  Wcc 
dXlc  Aic  dX2c 

=  A(l  +  e)  +  B- 


=  HAic  +  B  —  — 2 — 

^lc 

A  B 

(l  +  £)2  -  (l  +  e)3> 


B 

K 


(3.9) 


where  A  and  B  are  appropriate  constants  depending  on  the  values  of  C i, . . .  ,65.  Note  that 
in  this  calculation  we  use  the  relationship  A2C  =  A3C  =  y=  in  the  incompressible  CC  system 
under  tensile  deformation  in  the  X\  direction.  A  Taylor  expansion  of  the  negative  powers  of 
(1  +  e)  in  (3.9)  yields 


d  Wcc  X  2cdWa 


r*  .  I 


(3.10) 


Next  we  turn  to  the  contribution  of  Wpc  to  the  engineering  stress,  where  the  nonlinear 
assumption  (3.1)  modifies  the  results  of  the  linear  case.  We  find  that 


We  assume  that  £\  «  e  and  use  the  following  approximations 


/( 


A  ip 

Ale 


/(i±£i)=/(l)+/'(l)^  +  /"(1) 


1  +  £ 


l  +  £ 


£1  —  £ 
1  +  £ 


+  . 


/( 1)  +  /'W/rA  “  +  /'(1)(-£)(1  -  £  +  <E2  -  £S  +  .  .  .) 

1  +  £ 

OL 0  +  OL\£  +  OL2£"  +  .  .  .  , 


/(a 


Ale  n 

\J  Aip' 


=  /( 


1  +  £ 
1  +  £\ 
/"( i) 


)  «  /(Vl+i)  =  /(l)  +  /'(l)(vT+i -  1) 

(a/TTs  —  l)2  +  •  •  •  =  A)  +  (3\£  +  (32£2  +  /?3£T3  +  .  .  .  , 


where  we  expanded  a/1  +  £  in  powers  of  e.  With  similar  expansions  for  A|  ,  ^  A 3p2, 


iP 


Aip3^2,  Alc3^2  we  obtain  that 


A2C  ^/A2p,  i9Wpc 


Alp,  <9WpC  '*zc  r/''Z.p 

Aic  <9 A  ip  Aic  X2cJ  d\ 2p 


=  <5i£i  +  82e\  +  ...  +Ji£  +  72£2  +  •  •  • 


+h.o.t.(£k£{,  k,j  >  1).  (3-11) 

Combining  (3.10)  and  (3.11)  with  (3.8),  we  have  the  following  approximation  for  the  engi¬ 
neering  stress 

c \  (s,  £1)  —  C\£  -f-  c2£2  T  c3£3  T  h.o.t.  (sk ,  k  4)  -|-  hi£i  T  h.o.t.  (£^,  k  ^  2) 

+  h.0.t.(£fc£-(,  k,j  >  1)  ~  C]£  +  C2£2  +  C3£3  +  h]£i.  (3-12) 

We  remark  that  the  constitutive  equation  cq  =  ge(£)  +  Cx>£  +  /xi£i ,  where  ge  is  a  cubic 
polynomial  of  the  form 

QeijX)  =  Ci£  +  C2£2  +  C3£3,  (3.13) 


was  assumed  based  on  phenomenological  arguments  in  [10]  and,  as  we  shall  note  below  in 
Section  5,  it  has  been  successfully  used  in  reproducing  experimental  data  for  shear  deforma¬ 
tions  of  rubber  samples  as  well  as  for  large  tensile  deformations  of  a  rubber  rod.  We  also 
point  out  that  the  presence  of  the  e  term  in  this  formulation  for  o\  represents  an  approxi¬ 
mation  to  a  damping  mechanism  (certainly  present  in  the  case  of  highly  filled  rubber)  and 
does  not  result  from  the  nonlinear  reptation  formulation. 


4  General  tensile  models  with  multiple  relaxation  times 

The  above  model  with  (2.10),  (3.2)  (with  /  independent  of  t)  and  (3.12)  forms  the  basis  of 
a  general  class  of  nonlinear  deformation  models  for  polymers  in  tension.  More  generally,  the 
system  written  in  distributional  or  generalized  sense  (see  [7,  10])  has  the  form 

(4-i) 
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where 

Si(<S,  6",  <Si)  (J elastic i  ^visco^ l) 

(4.2) 

and 

(4.3) 

However,  such  models  are  based  on  the  tacit  assumption  that  all  molecules  relax  with  the 
same  relaxation  time  r.  There  is  substantial  experimental  evidence  [2,  16,  23,  24,  27]  to 
suggest  that  the  assumption  of  a  uniform  relaxation  time  is  not  valid.  Indeed,  when  fitting 
“elastic”  moduli 

E(u)  =  E'(u)  +  jE"(u)  =  f  ^  T  0(r)dr  +  j  j  0(r)dr 

Jo  1  +  C 0ZTZ  Jo  1  + 

to  response  data  in  the  frequency  domain,  it  is  often  necessary  to  use  a  probability  density 
(j)  for  the  distribution  of  relaxation  times. 

In  the  context  of  the  reptation  models  developed  here,  this  implies  that  one  should 
replace  the  engineering  stress  Si  in  (4.2)  by  a  probability  measure  or  probability  distribution 
dependent  stress-strain  law  of  the  form 

r  oo 

Si(M,  $)  =  creiast(£,£)  +  /  £i(r)d<f>(r),  (4.4) 

JTq 

where  t  — >  £\{t\  r)  is  the  solution  of  (4.3)  for  a  given  r,  <f>  is  a  probability  distribution  for  the 
relaxation  times  r,  and  To  >  0  is  a  lower  bound  on  possible  relaxation  times.  For  the  special 
case  of  an  absolutely  continuous  distribution  this,  of  course,  reduces  to 

/»oo 

Si(e,£,  <f>)  =  aelast(£,i)  +  /  £i (r)0(r)dr, 

J  TQ 

where  0  =  $'.  For  a  distribution  consisting  of  a  finite  number  of  Dirac  measures  with  atoms 
at  Ti,  72, . . . ,  tm,  respectively,  we  have 

M 

Si(£,d,<f>)  =  creiast(£,£)  +  J2pi£1(ri),  (4.5) 

1=1 

where  Pi  =  1  and  t  — >  £\{t\  Ti)  satisfies  (4.3)  with  r  =  r,,  i  =  1,  2, . . . ,  M. 

A  general  theory  of  existence  and  uniqueness  that  applies  to  systems  (4.1)  with  various 
approximations  of  (4.3),  and  (4.5),  i.e. ,  the  discrete  measure  case,  can  be  found  in  [1],  A 
theory  for  continuous  dependence  of  solutions  with  respect  to  parameters  (specifically  with 
respect  to  the  discrete  measure  in  (4.5))  can  be  pursued  in  the  context  of  the  Prohorov 
metric  topology  on  the  space  of  probability  distributions.  In  this  situation  this  is  equivalent 
to  the  weak-star  topology  on  the  distributions  if  they  are  viewed  as  a  subset  of  the  dual 
space  of  C[ro,r],  the  space  of  bounded  continuous  functions  on  a  finite  interval  [r0,r]  of 
relaxation  times  (see  [4,  5]  for  a  discussion  of  the  Prohorov  metric  and  its  use  in  inverse 
problems  for  systems  containing  uncertainty  in  the  “parameters”  in  the  setting  of  measure 
or  probability  distribution  dependent  systems).  Efforts  are  currently  underway  to  develop 
a  rigorous  theory  for  wcll-posedness  (including  continuous  dependence  with  respect  to  the 
measures  <f>)  of  systems  with  (4.4)  as  the  general  measure  dependent  engineering  stress  law. 
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5  Application  to  experimental  data 

Systems  such  as  (4.1),  (4.3)  with  (4.4)  certainly  pose  conceptual,  theoretical  and  computa¬ 
tional  challenges  even  when  treating  forward  or  simulation  problems.  Even  more  difficulties 
are  presented  by  estimation  or  inverse  problems  wherein  one  attempts  to  use  experimental 
data  to  determine  the  distribution  $,  parameters  in  aeiast,  or  /  in  (4.4)  and  (4.3),  respec¬ 
tively.  However,  significant  progress  on  such  problems  has  been  achieved.  In  a  recent  sum¬ 
mary  [7],  results  are  reported  on  using  E]  of  the  form  (4.5)  in  experimental,  computational 
and  theoretical  investigations  for  filled  rubber  rods  in  tensile  and  shear  deformations.  In 
the  following  we  describe  the  use  of  models  with  approximations  such  as  in  (3.3)  and  (3.12), 
that  is,  dropping  the  h.o.t.  terms,  to  fit  data  from  dynamic  experiments  for  rubber  samples 
in  tension. 

We  summarize  here  some  of  the  results  obtained  using  dynamic  experiments  with  a  rubber 
rod  with  a  tip  mass  at  one  end  in  uniaxial  tensile  deformation.  Similar  experimental  efforts 
to  validate  models  for  filled  rubber  samples  in  shear  were  successfully  carried  out  as  well  (see 
[7]).  The  experimental  device  constructed  specifically  for  these  tensile  validation  tests  at  the 
Thomas  Lord  Research  Center  of  the  Lord  Corporation  is  depicted  schematically  in  Figure 
2.  This  experiment  produced  data  {^}  consisting  of  time  measurements  of  force  (at  the  tip 


Figure  2:  Schematic  of  the  experimental  device  for  tensile  validation  tests 

of  the  rod,  x  =  0)  collected  by  a  load  cell.  The  data  corresponds  to  the  engineering  stress 
Ei  at  times  t,  at  the  top  of  the  rod,  x  =  0,  multiplied  by  the  cross  sectional  area  Ac.  Thus 
a  least  squares  formulation  for  estimation  of  parameters  q  has  the  form 

J(q )  =  5Z  I zi~  ASi(e(L,0),£(L:,0),<h)|2,  (5.1) 
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with  e(t,  x )  =  | ^{t,  x ),  where  Si  is  given  by  (4.5).  Here  £±  satisfies  (4.3)  with  /  approximated 
as  in  (3.3)  for  each  t*,  and  creiast  is  given  by 


(Telast(£,£ )  =  ge  (e)  +  CD£, 

where  ge  is  a  cubic  polynomial  as  in  (3.13).  To  be  more  precise,  experimental  data  suggested 
that  the  nonlinearity  /  was  not  the  same  when  tensile  deformation  was  increasing  as  when  de¬ 
creasing.  Hence  the  approximation  in  (3.3)  was  employed  with  two  sets  of  y/s,  one  set  {7 *nc} 
for  increasing  deformations  and  one  set  {y/ec}  for  decreasing  deformations.  The  parameter  q 
to  be  estimated  from  the  data  using  (5.1)  thus  consisted  of  p,  Cd,  {7i"c},  {7 fec},  <7,  c2,  C3  and 
the  Tj’s  in  (4.5).  For  the  results  described  here  we  fixed  M  —  2  in  (4.5)  with  p\  =  P2  =  1/2. 
For  highly  filled  rubber  we  found  estimation  with  one  uniform  relaxation  time  (M  =  1) 
would  not  adequately  describe  the  data.  The  details  for  the  case  M  —  2  are  given  in  [7,  10] 
and  a  typical  comparison  of  the  optimized  fit  model  (i.e.,  model  with  estimated  parameters) 
to  data  is  depicted  in  Figure  3. 


highly  filled  rubber  sample 


Figure  3:  Approximation  with  two  internal  variables  (dashed  line)  and  experimental  data 
(solid  line) 

Partial  theoretical  foundations  for  such  inverse  problems  are  available  [9]  with  a  complete 
theory  currently  being  pursued  by  the  authors  of  this  paper.  As  we  noted  above,  in  [1], 
wcll-posedness  results  (existence,  uniqueness)  for  forward  systems  of  the  form  (4.1),  with 
various  approximations  of  (4.3),  and  (4.5)  are  given  under  quite  general  assumptions  on  the 
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nonlinearities.  A  careful  formulation  of  the  associated  inverse  problem  for  estimation  of  <f> 
in  a  class  of  probability  measures  V  is  given  in  [5]  in  the  context  of  the  Prohorov  metric 
topology  on  V .  Computational  approaches  for  problems  similar  to  these  are  discussed  in  [4]. 
Our  efforts  to  develop  a  complete  theoretical  as  well  as  computational  framework  for  inverse 
problems  entailing  (4.1),  (4.3),  (4.4)  are  underway. 
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